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Abstract. A brief overview of selected topics in the theory and phenomenology of neutrino 
oscillations is given. These include: oscillations in vacuum and in matter; phenomenology 
of 3-fiavour neutrino oscillations and effective 2-flavour approximations; CP and T violation 
in neutrino oscillations in vacuum and in matter; matter effects on i/ M «-> v r oscillations; 
parametric resonance in neutrino oscillations inside the earth; oscillations below and above 
the MSW resonance; unsettled issues in the theory of neutrino oscillations. 



1. A bit of history... 

The idea of neutrino oscillations was first put forward by Pontecorvo in 1957 1 . Pontecorvo 
suggested the possibility of v <-> v oscillations, by analogy with K°K° oscillations (only one 
neutrino species - v e - was known at that time). Soon after the discovery of muon neutrino, 
Maki, Nakagawa and Sakata 2 suggested the possibility of neutrino flavour transitions (which 
they called "virtual transmutations"). 




Figure 1. Bruno Pontecorvo (1913 - 1993), Shoichi Sakata (1911 - 1970), Ziro Maki (1929 
2005) and Masami Nakagawa (1932 - 2001). 



2. Theory 

2.1. Neutrino oscillations in vacuum 

Neutrino oscillations are a manifestation of leptonic mixing. For massive neutrinos, weak 
(flavour) eigenstates do not in general coincide with mass eigenstates but are their linear 

1 Talk given at the XXII International Conference on Neutrino Physics and Astrophysics "Neutrino 2006" , Santa 
Fe, June 13-19, 2006 

2 On leave from the National Research Center "Kurchatov Institute" , Moscow, Russia 



combinations; the diagonalization of the leptonic mass matrices leads to the emergence of the 
leptonic mixing matrix in the expression for the charged current interactions, the relevant part 
of the leptonic Lagrangian being 

- C w+m = -^j={e. L ^ ll Vl\J L v L )W^ + diag. mass terms. (1) 



Here Vl and Ul are the left-handed unitary transformations that diagonalize the mass matrices 
of charged leptons and neutrinos. The leptonic mixing matrix U = V^Ul relates the left-handed 
components of the neutrino mass eigenstates and flavour eigenstates according to 

ka)=E C/ ail^ aSS ) (a = e,^r ] i = 1,2,3) . (2) 



For relativistic neutrinos, the oscillation probability in vacuum is 



P(v a -> u b ; L) 



2 



(3) 



For 2-flavour (2f ) oscillations, which are a good first approximation in many cases, one has 

\u e ) = cos 6 \v\) + sin6> \v2) , (4) 

\vn) = - sin# \vi) + cosO |z/ 2 ) , (5) 
and eq. (1) yields the 2f transition probability 

P tr = sin 2 2^ sin 2 Lj . (6) 

The modes of neutrinos oscillations depend on the character of neutrino mass terms: 

• Dirac mass terms (Pl^d^r + h.c): active - active oscillations v a L <-> VbL («> b = e, fx, r) 
Neutrinos are Dirac particles. 

• Majorana mass terms {pLfnL(yL) c + h.c): active - active oscillations v a i <-> u b L- 
Neutrinos are Majorana particles. 

• Dirac + Majorana mass terms (z^ltid-^r + ^l^l(^l) c + NrM(Nr) c + h.c): active - active 
oscillations v a L VbL\ active - sterile oscillations v a L (Nbn) c = (N£)l. 

Neutrinos are Majorana particles. 

Would an observation of active - sterile neutrino oscillations mean that neutrinos are Majorana 
particles? Not necessarily! In principle, one can have active - sterile oscillations with only Dirac 
- type mass terms at the expense of introducing additional species of sterile neutrinos with 
opposite lepton number L. 



2.2. Neutrino oscillations in matter ~ The MSW effect 

Matter can change the pattern of neutrino oscillations drastically. In particular, a resonance 
enhancement of oscillations and resonance flavour conversion become possible (Wolfenstein, 
1978; Mikheyev & Smirnov, 1985 |Sj). 

Matter effect on neutrino oscillations is due to the coherent forward scattering of neutrinos 
on the constituents of matter (fig. The neutral current interactions, mediated by the Z° 
boson exchange, are the same for active neutrinos of all three flavours (modulo tiny radiative 
corrections) and therefore do not affect neutrino oscillations. In contrast to this, charged current 



Figure 2. Lincoln Wolfenstein, Stanislav Mikheyev and Alexei Smirnov 




Figure 3. Neutrino 
interactions in matter 



interactions, mediated by the exchanges, are only possible for electron neutrinos because 
there are no muons or tauons in normal matter. This yields an effective potential of the electron 
neutrinos 

V e cc = V = V2G F N e , 

which leads to a modification of the nature of neutrino oscillations in matter. The 2f neutrino 
evolution equation in matter is 



d 
;— 
dt 



^fcos29 + V ^ sin 29 
~ cos 20 



A " r ' sin 20 A!: >" 



4E 
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(7) 



The mixing angle in matter 9 m , which diagonalizes the Hamiltonian on the r.h.s. of eq. l|7jl. is 
different from the vacuum mixing angle 9: 



sin 2 29 r . 



sin 2 29- i^ff 



cos 29 - V2G F N e } 2 + (%^) 2 sin 2 29 



The flavour eigenstates can now be written as 



cos C/ m \vi m 
— sin m I v\ 



If!/ 



2E 

+■ sin6» m \v 2m ) , 

+ cos e m \v 2m ) 



(8) 



(9) 
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where \v\m) and \v2m) are the eigenstates of the neutrino Hamiltonian in matter (matter 
eigenstates). The Mikheyev - Smirnov - Wolfenstein (MSW) resonance condition is 



V2G F N e 



Am 
2E 



cos 29 . 



(11) 



At the resonance u m — 45° (sin 2 26» m = 1), i.e. the mixing in matter becomes maximal. 

If the matter density changes slowly enough (adiabatically) along the neutrino trajectory, 
neutrinos can undergo a flavour conversion (see fig. 0J. In the adiabatic regime the transitions 
between the matter eigenstates \vi m ) and \v2m) are strongly suppressed, i.e. these states evolve 
independently. However, their flavour composition, which is determined by the mixing angle 9 m , 
varies with density: 9 m (N e » (iV e ) re s) ~ 90°, 9 m (N e = (N e ) res ) = 45°, 9 m (N e < (Ay rcs ) ~ 9. 
Therefore the state produced at high densities as, e.g., v e w v^m wm en d up at low densities as 




Figure 4. Adiabatic neutrino flavour 
conversion. Solid curves show the energy 
levels of neutrino matter eigenstates, 
dashed curves illustrate level crossing 
in the absence of mixing. Black and 
white filling corresponds to the weights 
of neutrino flavour eigenstates in given 
matter eigenstates. 
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a superposition of u e and with the weights sin 2 9 and cos 2 9, respectively. The adiabaticity 
(slow density change) condition can be written as 



sin 2 29 Am 2 
cos 29 ~2E 



-L p » 1 , 



(12) 



where L p - electron density scale height: L p = Kl/N^dNe/dxl^ 1 . 

A simple and useful formula for 2f conversion probability, averaged over production/detection 
positions (or small energy intervals), was derived in 0]: 



\ = ^ " ^ cos 2ft cos 20/(1 - 2P 1 ) 



(13) 



Here B% and 9f are the mixing angles in matter in the initial and final points of the neutrino path, 
and P' is the hopping probability, which takes into account possible deviations from adiabaticity: 
In the adiabatic regime P' <C 1, whereas in the extreme non-adiabatic regime P' = sin 2 (ft — 9f). 
The evolution equation for the neutrino system can be also written as 



where S = {Re(u*u^) , Im(v*v p ) , u*v e - 1/2} , 



(14) 



B = {(Am 2 /4E) sin20 TO , 0, V/2 - (Am 2 /4S) cos 29 m } . 



(15) 



The first equation here coincides with the equation for spin precession in a magnetic field. This 
analogy can be used for a graphical illustration of neutrino oscillations in matter (see fig. |SJ). 






Figure 5. Analogy between neutrino oscillations in matter and spin precession in a magnetic 
field. Left panel: oscillations in constant-density matter, middle and right panels - adiabatic 
and non-adiabatic conversions in matter of varying density (adopted from [3] ) . 



Another analogy of neutrino flavour conversion in matter is provided by a system of two 
coupled pendula |Bj (see fig. EJ). When the right pendulum gets a kick, it starts oscillating, 
but the left pendulum is almost at rest because the eigenfrequencies of the two pendula are 
very different. With the length 1% of the right pendulum slowly decreasing, its eigenfrequency 
approaches that of the left one, and when I2 = l\ the two frequencies coincide (the resonance 
occurs): both pendula oscillate with the same amplitude. When the length of the right pendulum 
decreases further, the amplitude of its oscillations decreases too, while the left pendulum starts 
oscillating with a large amplitude. This adiabatic transfer of the oscillation energy from one 
pendulum to another is analogous to the adiabatic neutrino flavour conversion. 
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Figure 6. Mechan- 
ical analogue of neu- 
trino flavour conver- 
sion in matter - two 
coupled pendula of 
variable lengths. 



Analysis of the solar neutrino data and the results of the KamLAND and CHOOZ reactor 
neutrino experiments has convincingly demonstrated that the (large mixing angle) MSW effect is 
responsible for the flavour conversion of solar neutrinos, thus resolving the long-standing problem 
of the deficiency of the observed flux of solar neutrinos. This is illustrated by the analysis of 
the Bari group, in which the strength of the matter-induced potential of electron neutrinos was 
considered a free parameter (fig. [7J). For more on MSW effect, see the talk of A. Friedland 8 . 



matter effects with standard size (V = aZ G F N e ) confirmed 



Bounds on a ual (Solor + CHOOZ + KamLAND) 




a MSW 



V(x)a a^Vtx) 



Figure 7. Results of the analysis 
of the solar, CHOOZ and KamLAND 
data with the standard matter-induced 
potential rescaled by a factor cimsw-, 
treated clS EL free parameter. The value 
omsw ~ 1 is strongly favoured [Zj. 



3. Phenomenology 

All the available neutrino data except those of the LSND experiment can be explained in terms 
of oscillations between the 3 known neutrino species - u e , j/„ and v T . If the LSND results 
are correct, they would most likely require the existence > 1 light sterile neutrinos v s (though 
some exotic scenarios also exist: CPT violation, violation of Lorentz invariance, mass- varying 
neutrinos, shortcuts in extra dimensions, decaying neutrinos, etc.). The MiniBooNE experiment 
was designed to confirm or refute the LSND claim, and the results are expected very soon. One 
should remember, however, that even if the LSND result is not confirmed, this would not rule 
out the existence of light sterile neutrinos and f a <-> u s oscillations, which is an intriguing 
possibility with important implications for particle physics, astrophysics and cosmology. From 



now on I will concentrate on 3-flavour (3f) oscillations of active neutrinos. For more on sterile 
neutrinos, see the talk of A. Kusenko [2]. 

3.1. 3-flavour neutrino mixing and oscillations 

For 3 neutrino species the mixing matrix depends in general on 3 mixing angles #12, $23 and 
#13, one Dirac-type CP-violating phase c>cp, and two Majorana-type CP violating phases o\^- 
The Majorana-type phases can be factored out in the mixing matrix according to U = UqK, 
K = diag(l , e tai , e t<72 ). The factor K does not affect neutrino oscillations and so will be omitted 
hereafter. Renaming Uq — > U, the relevant part of the leptonic mixing matrix can be written in 
the standard parameterization as 

(C12C13 S12C13 s 13 e~ t5cp \ 

-S12C23 - Cl2Si 3 S23e l,5cp C12C23 - Si 2 Si 3 S23e l,5cp C13S23 , (16) 

S12S23 - Ci2Si 3 C2 3 e l<5cp -C12S23 - Sl2Sl3C23e l<5cP C13C23 / 

where Sy = sin 0y and = cos 8ij . 

Neutrino oscillations probe the neutrino mass squared differences, which satisfy Am s 2 ol = 
Am^ <C Am| 2 — Ami^ = Am^ tm . Accordingly, there are two possible orderings of the neutrino 
masses: normal hierarchy, when the mass eigenstate 1/3, separated from v\ and v<i by the largest 
mass gap, is the heaviest one, and inverted hierarchy, when v% is the lightest state (see fig. |S|). 



Normal hierarchy: 



Inverted hierarchy: 





Figure 8. Normal 
and inverted neu- 
trino mass order- 
ings. The different 
fillings show the rel- 
ative weights of dif- 
ferent flavour eigen- 
states in given mass 
eigenstates. 



3. 2. 2f and effective 2f approximations 

In many cases 2f description of neutrino oscillations gives a good first approximation. The 
reasons for this are (i) the hierarchy of Am 2 : Am^ <C Am^ tm , and (ii) the smallness of \U e s\. 
There are exceptions, however: when oscillations due to the solar frequency (oc Am 2 ol ) are not 
frozen, the probabilities P(y^ v T ), P{ u n v a) an d P(y T — ► v T ) do not have a 2f form 
jlOj . However, even for the probabilities of oscillations involving z/ e , the corrections due to 3- 
flavourness can be as large as ~ 10%, i.e. are at the same level as the accuracy of the present-day 
data, and so cannot be ignored. In addition, there is a number of very interesting pure 3f effects 
in neutrino oscillations. Therefore, 3f analyses are now a must. 

For oscillations driven by AmL, the third neutrino mass eigenstate 1/3 essentially decouples. 
However, there is still a "memory" of this state through unitarity, which results in the emergence 
of various powers of C13 in the expressions for transition probabilities. An example is the survival 
probability of solar v e 11 (the same expression also applies for the survival probability of reactor 
v e observed in KamLAND) : 



P D (Ve -> Ve) - cf 3 P 2ee (Am2 1 , #12, C13V) + sf 3 



(17) 



The term sf 3 is tiny and can be safely neglected. Another example is the day-night effect for 
solar v e : while the day-time survival probability Pryiye) oc cf%, the difference of the night-time 
and day-time probabilities Pjssiy^) — Po(V e ) oc cf 3 ^2jE]- Deviations from 2f expressions (which 
correspond to the limit #13 —* 0) may be substantial: for the maximal currently experimentally 
allowed values of sin 2 26*13 one has (1 — cf 3 ) ~ 0.1, (1 — cf 3 ) ~ 0.13. 

For oscillations of reactor v e) the survival probability can to a very good accuracy be written 

as 

P- e - e ~ 1 - sin 2 26 13 ■ sin 2 - 4 ^ 29 12 • sin 2 . (18) 

Since the average energies of reactor antineutrinos E ~ 4 MeV, for experiments with relatively 
short baseline (L < 1 km), such as CHOOZ, Palo Verde and Double CHOOZ, one has 
(Am^/AE) L<1. Eq. (HE) then reduces to 

P(y e —>■ i> e ; L) = 1 — sin 2 20 13 • sin 2 f , (19) 

i.e. takes the 2f form. Note that the "solar" term ~ sin 2 29\2 in ()18|) cannot be neglected if 
#13 ^ 0.03, which is about the reach of currently discussed future reactor experiments. 

For the unique long-baseline reactor neutrino experiment KamLAND (L ~ 170 km) one has 
(Am'^i/'iE) L > 1, (Am 31 /4P) L>1, and the v e survival probability takes the effective 2f form 
as in eq. (|17j) . Note that matter effects in KamLAND should be of order a few per cent, i.e. can 
be comparable with 3f corrections due to #13 7^ 0. 

3.3. Genuine 3f effects in neutrino oscillations 

These are, first of all, CP and T violation. CP violation results in P{y a — ¥ ^b) 7^ P(v a ^>)> 
whereas T violation leads to P{y a — > u^) ^ P{yb — > v a ). Under the standard assumptions of 
locality and normal relation between spin and statistics, quantum field theory conserves CPT. 
CPT invariance of neutrino oscillations in vacuum gives P{y a — > v^) = P{pb ~^ v a ). Therefore 
CP violation implies T violation and vice versa. 

One can consider the following probability differences as measures of CP and T violation: 

AP a c b P = P{ya ^ Vb ) _ P{ p a ^ 9b) ; AP t = P{ya Vh) _ P{yh Va) _ (20) 

From CPT invariance, for oscillations in vacuum one has 

AP a c b P = A pT ) AP fl c a P = . (21) 

In the 3f case there is only one Dirac-type CP-violating phase <5cp and therefore only one CP 
and T violating probability difference: 

AP e C p = APCP = AP T c e p = AP, where (22) 



AP = - 4si2 C12 S13 c 2 3 S23 c 23 sin 5 C p 

, Am 12 A • f Am h r\ ■ /Am 2 _ 

sm{ ^f L +sm h# L +sm h# L 



(23) 



This probability difference vanishes when one or more of the following conditions are satisfied: 
at least one Amfj = 0; at least one #y = or 90°; 5qp = or 180°; in the regime of complete 
averaging; in the limit L —* (AP — > as L 3 ). Obviously, the effects of CP and T violating are 
very difficult to observe! For more on that, see the talk of O. Mena |14j . 



CP violation and T violation in v oscillations in matter. Normal matter (with number of 
particles 7^ number of antiparticles) violates C, CP and CPT, which leads to a fake (extrinsic) 
CP violation in neutrino oscillations. It exists even in the 2f limit and may complicate the study 
of the fundamental (intrinsic) CP violation. 

The situation with T- violation in matter is different: matter with density profile symmetric 
w.r.t. the midpoint of neutrino trajectory does not induce any fake T violation. Asymmetric 
profiles do, but only for N > 2 flavors |15l I16j . Matter-induced T violation is an interesting 
pure 3f effect; it may fake fundamental T violation and complicate its study (extraction of Sep 
from the experiment). However, it is absent when either U e 3 = or Am? ol = (2f limits) and 
thus is doubly suppressed by both these small parameters. Therefore its effects in terrestrial 
experiments are expected to be very small |16j . 

Matter effects on ^ v T oscillations. In the 2f limit, matter does not affect «-> u T 
oscillations (because the matter-induced potentials V{u^) and V(u T ) coincide up to tiny radiative 
corrections). However, this is not true in the full 3f framework 17 . In particular, for oscillations 
inside the earth there are ranges of baselines and neutrino energies for which the matter effect 
can be very large (fig. 03 left panel, E ~ 5 — 10 GeV). If one ignores them, one may end up 
with a negative expected flux of oscillated z/„ in atmospheric neutrino experiments (fig. 03 right 
panel) . 

Figure 9. Left panel: P^. 
Right panel: oscillated flux 
of atmospheric v^. Am^ = 
2.5 x 10~ 3 eV 2 , sin 2 #13 = 
0.026, 6*23 = vr/4, Ara^ = 0, 
L = 9400 km. Red (dark) 
curves - with matter effects, 
green (light) curves - without 
matter effects on P UT . 



3-4- Parametric resonance in neutrino oscillations in matter 

The MSW effect is not the only possible way matter can influence neutrino oscillations. Another 
interesting possibility is a parametric enhancement of neutrino oscillations in matter |18[ I19j . 
Parametric resonance in oscillating systems with varying parameters occurs when the rate of the 
parameter change is correlated in a certain way with the values of the parameters themselves. A 
well-known mechanical example is a pendulum with vertically oscillating suspension point (fig. 
I10p : when the frequency O and amplitude A of these oscillations are in a special correlation with 
the eigenfrequency u and amplitude a of the pendulum, the pendulum can turn upside down and 
start oscillating around the vertical, normally unstable, equilibrium point. Neutrino oscillations 






Figure 10. Parametric resonance in oscilla- 
tions of a pendulum with vertically oscillating 
point of support. For small-amplitude oscilla- 
tions the resonance condition is f2 res = 2u)/n 
(n = 1.2.3...). 



in matter can undergo parametric enhancement if the length and size of the density modulation is 
correlated in a certain way with neutrino parameters. This enhancement is completely different 
from the MSW effect; in particular no level crossing is required. An example admitting an 
exact analytic solution is the "castle wall" density profile |19| 120] (see fig. Illj) . The resonance 
condition in this case can be written as [20] 

A3 = —(sin (fti cos 02 cos 26i m + cos 4>\ sin cos 202m) = , (24) 

where cpi^ are the oscillation phases acquired in layers 1 and 2 and 6 m i,2 are the corresponding 
mixing angles in matter. 




Figure 11. Parametric resonance in the case of a "castle wall" density profile. Coordinate 
dependence of the potential V (left panel) and of the transition probability P (right panel). 

The earth's density profile seen by neutrinos with core-crossing trajectories can be well 
approximated by a piece of this castle wall profile (fig. IT2*|) . Interestingly, the parametric 




Figure 12. The earth density profile according 
to the PREM model [2^ and its approximation by 
a piece of the "castle wall" profile. 



Lpcm] 

resonance condition (|24|) can be satisfied for oscillations of core-crossing neutrinos in the earth 
for a rather wide range of zenith angles both at intermediate energies |221 I2.H1 I2()| and high 
energies |24| (see figs. II HI I14|) . The parametric resonance of neutrino oscillations in the earth 
can be observed in future atmospheric or accelerator experiments if #13 is not too much below 
its current upper limit. 

3.5. Some recent developments 

When V <C Am 2 /2E (oscillations of low-E neutrinos in matter or, equivalently, oscillations 
in low-density matter), matter effects on neutrino oscillations are small and can be considered 
in perturbation theory. This gives simple and transparent formulas describing, in particular, 
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Figure 13. Parametric 
resonance in oscillations 
of core-crossing neutri- 
nos. Left panel: cos = 
-0.89, sin 2 20i3 = 0.01; 
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0.15, small window shows 
the values of A3. 
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Figure 14. Coordinate depen- 
dence of the transition probability 
P2 = 1 — P ee (upper panels) and of 
the total adiabatic oscillation phase 
$ (lower panels) for a number of 
values of E and cos 0^. The shaded 
areas show the earth's density pro- 
files seen by the neutrinos. Figure 
courtesy of M. Maltoni. 



oscillations of solar and supernova neutrinos in the earth. The earth matter effects can be 
expressed through the regeneration factor f reg = P® e — P2e C > wriere ?2e is the probability for v 2 
to become v e upon traversing the earth. In the 3f framework one has |13j 



2e 



rjvac 
r 2e 



cf 3 sin 2 29 



12 



dx V{x) sin 



2 u(x')dx' 



(25) 



where 



lo(x) 



[cos 20i2 5 - cj 3 V(x) /2] 2 + 5 2 sin 2 26 12 , 6 



Am|i 
AE 



(26) 



The 2f (0i3 = 0) version of these equations was derived in [25] (see also [HI]). 

The regeneration factor as the function of the cosine of the nadir angle 2 is shown in 
fig. 1151 As can be seen in the left panel, in the case of perfect energy resolution one could 
expect a significant increase of the regeneration factor for core-crossing trajectories. However, 
experimentally no such an increase was observed in the cos Z dependence of the day-night 
signal difference for solar neutrinos, which is rather flat. As was shown in 25.}, this comes about 
because of the finite energy resolution AE of the detectors, which leads to a suppression of the 
effects of the earth density variations that are far from the detector (see the right panel of fig. 
IT3|) . the attenuation length being d ~ l osc (E/AE). 

Oscillations of high energy neutrinos in matter or, equivalently, oscillations in dense matter 
(V > 5 = Am 2 /AE), can also be very accurately described analytically. The transition 
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Figure 15. Left panel: earth regeneration factor for E = 10 MeV neutrinos, perfect energy 
resolution. Black (solid) curve - numerical calculation, red (dashed) curve - analytic result. 
Right panel: regeneration factors averaged over three different intervals of energy |26| . 



probability for oscillations in a matter of an arbitrary density profile is given by [2 



P = 5 2 sin 2 29 



L 2 



o 



4>{x) = dx u)(x ) . 
Jo 



(27) 



The accuracy of this approximation is quickly increasing with neutrino energy (see the right 
panel of fig. I13| where the exact results are shown by solid curves and the analytic results, 
by dashed curves). Eq. (|27j) also allows a simple analytic interpretation of the two prominent 
parametric peaks in the core region, seen in this figure 24 . 



4. Unsettled issues? 

The theory of neutrino oscillations is quite mature and well developed now. However, it is far 
from being complete or finished, and a number of basic questions are still being debated. Below 
I list some of these questions (given in italics), along with my short answers to them: 3 

• Equal energies or equal momenta? 

- Neither equal E nor equal p assumptions normally used in the derivations of the oscillation 
probability are exact. But for relativistic neutrinos, both give the correct answer. 

• Evolution in space or in time? 

- This is related to the previous question. For relativistic neutrinos both are correct and 
equivalent. Fortunately, in practice we only deal with relativistic neutrinos. In the non- 
relativistic case the very notion of the oscillation probability is ill-defined (the probability 
depends on both the production and detection processes). 

• Claim: evolution in time is never observed 

- Incorrect. Examples: K2K, MINOS (and now also CNGS) experiments, which use the 
neutrino time of flight in order to suppress the background. 

• Is wave packet description necessary? 

- Yes, if one wants to rigorously justify the standard oscillation probability formula. Once 
this is done, the wave packets can be forgotten unless the issues of coherence become 
important. 

• Do charged leptons oscillate? 

- No, they don't. 



3 Detailed discussion could not be given for the lack of time. 



• Is the standard oscillation formula correct? 

- Yes, it is. In particular, there is no extra factor of two in the oscillation phase, which 
is sometimes claimed to be there. However, it would be theoretically interesting and 
important to study the limits of applicability of the standard approach. 

A number of subtle issues of the neutrino oscillation theory still remain unsettled (e.g., 
rigorous wave packet treatment, oscillations of non-relativistic neutrinos, etc). At present, this 
is (rightfully) of little concern for practitioners. 

What are interesting future tasks for the theory and phenomenology of neutrino oscillations? 
These include the search for the best strategies for measuring neutrino parameters, study of 
subleading effects and effects of non-standard neutrino interactions and of the domains of 
applicability and limitations of the current theoretical framework. Future experimental results 
may also bring some new surprises and pose more challenging problems for the theory! 
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